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Abstract. We calculate the leading isospin conserving few-nucleon contributions to pion scattering on ^H, 
■^He, and *He. We demonstrate that the strong contributions to the pion-nucleus scattering lengths can be 
controlled theoretically to an accuracy of a few percent for isoscalar nuclei and of 10% for isovector nuclei. 
In particular, we find the yr-^He scattering length to be (62 ± 4 ± 7) x 10~^m^^ where the uncertainties 
are due to ambiguities in the vr-N scattering lengths and few-nucleon effects, respectively. To establish this 
accuracy we need to identify a suitable power counting for pion-nucleus scattering. For this purpose we 
study the dependence of the two-nucleon contributions to the scattering length on the binding energy of ^H. 
Furthermore, we investigate the relative size of the leading two-, three-, and four-nucleon contributions. 
For the numerical evaluation of the pertinent integrals, a Monte Carlo method suitable for momentum 
space is devised. Our results show that in general the power counting suggested by Weinberg is capable to 
properly predict the relative importance of A'^-nucleon operators, however, it fails to capture the relative 
strength of N- and (A'^ + l)-nucleon operators, where we find a suppression by a factor of 5 compared to 
the predicted factor of 50. The relevance for the extraction of the isoscalar vr-N scattering length from 
pionic and ^He is discussed. As a side result, we show that the calculation of the vr-^H scattering length 
is already beyond the range of applicability of heavy pion effective field theory. 
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1 Introduction 

The probably most fundamental quantities characterizing 
pion scattering off nucleons are the pion-nucleon (tt-N) 
scattering lengths. Obviously being low energy observ- 
ables, one should be able to understand their size based 
on the effective field theory (EFT) of Quantum Chro- 
■ modynamics (QCD), namely Chiral Perturbation Theory 
(ChPT). In fact, the impact of symmetries on these quan- 
tities is known for a long time FTIEI (for a recent review see 
Ref. [3]), but for a quantitative analysis of high accuracy 
it is of utmost importance to understand now in detail the 
higher order corrections, either isospin symmetric ones or 
isospin symmetry violating (IV) ones, based on ChPT, for 
they point directly at QCD dynamics. 

Despite an on-going effort to learn more on the scat- 
tering lengths experimentally, they are still not very accu- 
rately known [T. This is especially true for the isoscalar 
scattering length which is so small, that a high control 
over various other, normally neglected effects, like IV, is 
necessary for a reliable extraction of a^"*"-* from experi- 
ment. On the theoretical side important progress for IV 
TT-N scattering is reported in Refs. [51E1I7]. A first extrac- 



tion of the leading charge symmetry breaking (a special 
case of IV interactions) tt-N scattering amplitude is re- 
ported in Refs. [8lf9llT0] based on an analysis of the re- 
cent observation of the IV forward-backward asymmetry 
in pn — >■ c?7r° [llj . 

The most promising way to get experimental access to 
the TT-N scattering lengths is the investigation of pionic 
atoms, the simplest system being pionic hydrogen. Many 
efforts have been undertaken to measure energy shift and 
width of this system, allowing for an extraction of the 
isoscalar and isovector tt-N scattering lengths a^+^ and 
a^"^ assuming that isospin symmetry is exact up to stan- 
dard Coulomb interactions P^^ITB] . 

In principle, it appears to be possible to extract both 
the isoscalar and the isovector scattering length from pio- 
nic hydrogen data, for the measurement provides two in- 
dependent observables, namely the level shift compared to 
the QED value induced by the strong interaction as well as 
the width, where the transition iT~p — > 7r°n gives a signif- 
icant, known contribution. However, in order to improve 
the extraction accuracy and to better control systematics, 
an additional source of a^^^ is desireable. In this context, 
it is important to study more complex systems, like pionic 
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deuterium, '^He, or "^He, for the wave functions for those 
nuclei can be calculated with high accuracy. It is easy to 
see that the pion-nucleus scattering length Ut^-a can be 
written as 

-|-few-nucleon corrections 
+IV corrections (1) 

Here, A denotes the mass number of the nucleus, Tj, the 
third component of its isospin, and Q^r the mass and 
charge of the pion and ttzat the nucleon mass. We will not 
further discuss the IV corrections here, instead we refer to 
Refs. [Hill] for a review on the current status of IV in pion- 
nucleus scattering. Especially, the isospin T = nuclei 
and ^He yield important new information since contribu- 
tions from a^~^ are suppressed and the contribution of a^+^ 
is more prominent. From Eq. ([1} it becomes clear that an 
understanding of these more complex systems requires to 
understand with a high accuracy the few-nucleon contri- 
butions to the scattering length. For tt-^H scattering they 
have been studied for decades within potential multiple 
scattering theory, see e.g. Refs. jl5ul6ul7i and references 
therein. A first attempt to calculate few-body effects sys- 
tematically and model-independently, i.e. within ChPT, 
was done by Weinberg [TH]. He proposed to organize the 
transition operator as an expansion in mT^j A with A be- 
ing the chiral symmetry breaking scale of the order of 
1 GeV. This transition operator was then convoluted in 
Ref. [IH] with phenomenological models for the nucleon- 
nucleon (NN) interaction. This work has been refined in 
Ref. [l^. The explicit calculations showed that part of 
the few-nucleon corrections seem to be suppressed. Maybe 
even more troublesome for the extraction of tt-N scattering 
lengths is the observation that the few-nucleon corrections 
were somewhat larger than naive dimensional analysis pre- 
dicted. Such issues were addressed in [3D]. Here a refined 
power counting, the so-called (5~counting, was proposed 
that takes the small binding energy of ^H explicitly into 
account. Assuming that typical momenta of the nucleons 
are given by the binding momentum 7 = y^dm/v, with 
Es, denoting the deuteron binding energy, it was found 
that the estimated magnitudes of various few-body cor- 
rections changed drastically and thus the series needed to 
be reorganized compared to the original power counting 
used in |18j . Obviously, it is important to investigate the 
validity of this assumption, before extending the approach 
to even more complex systems or towards higher orders. 

As mentioned above, the second issue of the first calcu- 
lations of two- nucleon corrections to tt-^H scattering is the 
relatively large size of the leading few-nucleon terms. Re- 
cently, higher order contributions to this class of diagrams 
have been investigated in Ref. [30] and Refs. [3TJ[32]. Es- 
pecially, the results of the latter references indicate that 
the small scale 7 does not play a role in the higher order 
terms studied in Refs. [3T1I22) . It is however not clear at 
this point, whether the hierarchy of two-, three- and four- 
nucleon contributions is as expected by the power count- 
ing. An understanding of the systematics of such higher 



order corrections is mandatory when analyzing data on 
^He and "^He. 

In this paper, we study the presumably leading two-, 
three- and four-nucleon contributions explicitly. In Sec. [3J 
we will discuss the implications of using the different power 
countings in more detail. In Sec. [31 we repeat the well- 
known expressions for two-nucleon operators, derive a com- 
plete set of leading three-nucleon operators and give ex- 
pressions for the probably most relevant example of a 
four-nucleon operator. The contribution to the scattering 
length based on these expressions requires the evaluation 
of high dimensional integrals involving the operators and 
few-nucleon wave functions. Such a problem will also ap- 
pear in other applications. Therefore, we explain our nu- 
merical method for such calculations in detail in Sec. [H 
In Sec. ISl we turn to ^H and study the dependence of the 
leading two-nucleon contributions on the binding energy 
based on an unphysically over- or underbound ^H. The re- 
sults will shed light on the validity of Q-counting. We then 
introduce the input into our calculations for the required 
few-nucleon wave functions in more detail in Sec. [6l This 
is a prerequisite for our further studies. Sec. [71 is devoted 
to '^He where the scattering length of 7r-^He system is cal- 
culated. For the first time, we give explicit results for the 
contribution of three-nucleon corrections to pion-nucleus 
scattering. Finally, in Sec.jSl we estimate four-nucleon con- 
tributions to TT-'^He scattering based on an explicit calcu- 
lation. We then summarize our findings and conclude in 
Sec.jSl More technical aspects of this work are deferred to 
the appendices. 

2 Power-counting schemes for pion-nucleus 
scattering 

The central goal of this work is to provide the calcula- 
tions necessary to extract the isoscalar scattering length 
from data on the pion-nucleus scattering length. The basis 
for this is a reliable power counting which allows one to 
identify the relevant operators and provides a framework 
to estimate the diagrams not included. To reach maximal 
predictive power we include in each channel all diagrams 
that contribute up to one order lower than the contribu- 
tion of the leading, unknown counter term. Especially the 
estimate of higher order operators requires some care. We 
here estimate the contributions from higher order opera- 
tors, not calculated explicitly, by two methods: 

• We demonstrate that all A^-nucleon diagrams included 
can be organized in a series in powers of some param- 
eter X- Thus, a calculation performed up to order x" 
is expected to have an uncertainty of order x^""*"^^- 

• Using a set of wave functions generated for regulators 
varying over a large range of values, we investigate 
quantitatively the regulator dependence of the various 
matrix elements. This regulator dependence should be 
absorbed into a pertinent counter term. Therefore, we 
get access to the possible counter term contribution. 

As we will see, for isoscalar nuclei both methods give con- 
sistent results supporting our claim that the scattering of 
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pions off isoscalar targets at threshold can be calculated 
with an uncertainty of the order of 5 % of the leading 
two-nuclcon contributions: 1 x 10^"^ m^^ for tt-^H and 
2 X 10""^ m^^ for ir-^He scattering. For isovector nuclei 
the contribution of the leading contact term is estimated 
via scaling of the leading isoscalar contribution. This al- 
lows us to estimate the accuracy of the calculation for 
TT-'^He scattering to be of the order of 10%. 

It should also be stressed that up to now we are not 
able to identify a power counting that allows one to relate 
operators with different numbers of active nucleons, say 
one-nucleon with two- or three-nucleon operators. How- 
ever, although we do not quantitatively understand the 
reasons for successive suppression of operators with an 
increasing number of active nucleons, the suppression ob- 
served numerically is quantitatively sufficiently large to al- 
low for controlled results. In addition, it turns out that the 
successive suppression observed between N- and {N -1-1)- 
nucleon operators is the same for = 1, 2, 3, which is 
important for controlled uncertainty estimates. We also 
demonstrate that it is indeed possible to estimate reliably 
the relative importance of various 7V-nucleon operators. 

What we need to carry out this program is a reliable 
power counting. There are two scales relevant for pion- 
nucleus scattering at threshold, namely the pion mass, 
TOtt, and the nucleus binding momentum, 7, which may 
either be evaluated via 

7 = > (2) 

with n for the reduced mass of a single nucleon with re- 
spect to the remainder and e for the binding energy with 
respect to the first break-up channel, or via 

7 = y/2mN{E/A) , (3) 

with (E/A) for the binding energy per nucleon. For the 
nuclei of relevance here — ^H, ■^He, and ''He — both for- 
mulas give similar answers, namely 46 and 46 MeV, 82 and 
69 MeV, and 167 and 115 MeV for 7 and 7 for the three 
nuclei in order. In the following we will therefore only use 
the quantity 7 for the binding momentum. 

Let q denote some generic momentum. There are two 
different power-counting schemes proposed for pion-nucleus 
scattering at threshold. These differ in how q is treated rel- 
ative to niTT and 7. In Ref. [TB], referred to as Weinberg 
counting, the assignment 

q'^m^ (4) 

is proposed while in Ref. '20^, in the following called Q~ 
counting, 

q r^j <^mT, (5) 

is used. In [20] the authors argue that the latter assign- 
ment is to be preferred, since, for tt-^H scattering, i) it 
naturally explains why the contribution of Diagram (a) 
of Fig. [T] is more than an order of magnitude larger than 
those of (5) and (c); ii) it quantitatively explains the mag- 
nitude of the contribution of the diagram of Fig. [5] and 
in) it quantitatively explains the contribution of the boost 



correction — this term emerges as a recoil correction to 
the leading isoscalar tt-N scattering amplitude. 

As in the ttN sector, chiral symmetry enforces that op- 
erators that are odd under the exchange of the external 
pions — isovector operators — appear with at least one 
derivative acting on the pion field, while those that are 
even — isoscalar — have at least two derivatives or one 
insertion of the quark mass matrix. As a consequence, 
isoscalar operators are one order in niT^/mpf suppressed 
compared to their isovector counter parts. This is another 
argument in favour of using isoscalar targets to get high 
accuracy information on the tt-N scattering lengths: since 
the contribution of the leading order counter term sets the 
limit for the accuracy for the extraction of the tt-N scat- 
tering lengths, the analysis of the scattering off isovector 
nuclei provides values for the scattering lengths about an 
order of magnitude less accurate. 

The pertinent few-body corrections can now be classi- 
fied in terms of q and m^. Because of the arguments above, 
the leading isoscalar 4N7r7r counter term contributes to 
the TT-^H scattering length at order (mTr/m^)'^ , where for 
simplicity we identified the chiral symmetry breaking scale 
with the nucleon mass, while the leading two- nucleon cor- 
rection — Diagram (a) of Fig. [1] — scales as (rriT^/q)^. 
Explicit calculation gives that the leading two-nucleon di- 
agram contributes about 20x10~^to~^ to the tt-^H scat- 
tering length. Thus in the Weinberg scheme we estimate 
about Ix 10~^ for the leading counter term contribu- 
tion while in Q-counting the leading counter term is ex- 
pected to contribute only O.lx lO^'^ '^^^j since the leading 
two-nucleon term benefits from an enhancement of order 
(tott/?)^ — (™7r/7)^ in Q-counting. Thus again, if we aim 
at a reliable estimate of the accuracy of the calculations 
performed we need to understand which one of the power- 
counting schemes describes the hierarchy of diagrams for 
pion-nucleus scattering at threshold. This will be done in 
Sec. 01 

It turns out that Weinberg and Q-counting predict 
a very different binding energy dependence of ratios of 
few-body corrections. These relations can be tested em- 
pirically: the NN potential at leading order of the chiral 
expansion comprises l7r-exchange and two counter terms. 
One of them can be adjusted such that any given (physi- 
cal or unphysical) deuteron binding energy can be repro- 
duced. This has been done to investigate the binding en- 
ergy dependence of the relevant ratios of few-nucleon con- 
tributions. The main conclusion of Sec. l5.2l will be that the 
binding energy dependence of these ratios is in accordance 
with the prediction from Weinberg — up to logarithmic 
corrections — and in strong disagreement to the predic- 
tion from Q-counting. 

In addition to providing some explanation for the ap- 
parent suppression of a^^'^'^ vs. a^^°^, in Ref. [3D] it was 
argued that boosted tt-N amplitudes give a significant con- 
tribution (number in) in the list given above), in discrep- 
ancy to what is expected from Weinberg counting. On the 
other hand, since the corresponding operator is propor- 
tional to the square of the nucleon momentum, this would 
imply an observable effect of the nucleon kinetic energy in- 



4 



S. Liebig et al.: Towards a high precision calculation for the pion-nucleus scattering lengths 



side the deuteron, which is in conflict with general prop- 
erties of field theories. In Ref. [22] both problems were 
solved: once the Z\-isobar is included as an explicit degree 
of freedom, simultaneously the residual boost contribu- 
tion becomes sufficiently small to be in line with Wein- 
berg counting and appears at an order where there is also 
a counter term. 

Finally, we analyze the contribution of the triple scat- 
tering diagram of Fig. [2] in Sec. 15.31 In Weinberg count- 
ing this diagram is estimated to contribute only at or- 
der (tott/^jv)^ relative to the leading two-nucleon oper- 
ator. On the other hand numerical studies revealed that 
its actual value exceeds this estimate by about an order 
of magnitude. In Sec. 15.31 we will demonstrate that this 
enhancement is not due to the smallness of the deuteron 
binding energy, contrary to the claim in Ref. [20 , but due 
to an integrable singularity in the corresponding expres- 
sion that produces an enhancement by a factor of tt^. In 
view of this, also item ii) of the list that was used to argue 
in favour of Q-counting does not apply anymore. 

We are led to conclude that Weinberg counting gives 
a more reliable power counting scheme than Q-counting. 

With these remarks, we now want to look at the pos- 
sible few-nucleon contributions in more detail. 



3 Few-nucleon contributions 

Our main concern here are the few-nucleon corrections. 
For an evaluation of such contributions, expectation values 
of the basic pion-few-nucleon (tt-A) amplitudes with the 
few-nucleon wave functions need to be calculated. Such 
wave functions can be obtained using standard methods 
to solve the few-nucleon Schrodinger equation based on 
nucleon-nucleon (NN) interactions 23 . Whereas in the 
first attempts to understand such contributions, wave func- 
tions based on phenomenological NN interactions were 
used P^ITO] . it is by now standard to also employ wave 
functions generated with nuclear interactions based on 
ChPT. The basic idea is to use naive dimensional analysis 
for a potential, which is then used to solve a Schrodinger 
equation [^[^ . Several groups have developed NN inter- 
actions based on this approach pM^I^ and also three- 
(3N) and four-nucleon (4N) interactions have been formu- 
lated [29ll30] and employed [31], see also Ref. [S^ for a 
recent review. Obviously, employing chiral interactions is 
preferable to calculate the pion-nucleus scattering length 
since nuclear interactions and tt-A amplitudes will be con- 
sistent. 

In this work, we will present results on both kinds of 
interactions. We still show results based on phenomeno- 
logical interactions, namely AV18 [33j, Nijmegen 93 [34. . 
and CD-Bonn [3S]. For ^He and ^He, we augment the nu- 
clear Hamiltonian by a 3N interaction, so that the bind- 
ing energies of these nuclei are reasonably well reproduced 
[36] . These results may serve as benchmark and might give 
indications on the size of the model-dependence of older 
calculations. We also employ wave functions based on chi- 
ral nuclear interactions between the leading order (LO, or- 



der in the chiral expansion) and next-to-next-to-leading 
order (N^LO, order 3 in the chiral expansion). 

Chiral interactions require a regularization scheme in 
order to obtain a well-defined Schrodinger equation. Most 
realizations use a momentum cutoff of the order of 500 MeV 
to this aim. For the leading order, at least when restricting 
oneself to S-wave interactions, it is possible to obtain fits 
for a much larger range of cutoffs [37 . These attempts 
have triggered a controversy in the community 38 . It 
should be noted that an inconsistency becomes apparent 
for large cutoffs, when higher order NN interactions are 
used [22 ■ Here, we use a wide range of cutoffs only for 
LO interactions, so that also this inconsistency does not 
apply. This allows us to use results for a wide range of 
cutoffs to estimate the size of leading counter terms. Pa- 
rameter sets for the S'-wave contact interactions are given 
in Tables [31 S] and |SJ Since these calculations are by no 
means high precision ones, we have neglected the minor 
contribution of higher partial waves in this case. 

For the higher order chiral forces, we have employed 
order 2 (NLO) and order 3 (N^LO) ones of Ref. [28]. In 
the N^LO case, we added, as required by power counting, 
3N forces, which were tuned to reproduce the '^He binding 
energy and N-^H scattering lengths. 

Since we restrict ourselves to leading, tree-level n-A op- 
erators, we may derive them using Feynman diagrams as 
done below. The reduction from four-dimensional quan- 
tities to three-dimensional ones can be easily performed 
using the on-shell energies of nucleons and pions. The per- 
tinent integrals involving the wave functions and operators 
will be given below since their form depends on the num- 
ber of nucleons involved. 



3.1 Leading two-nucleon contributions 

The leading two-nucleon contributions are known for many 
years [l8lfT9l[20| . In the following, we call the numeri- 
cally most important contribution, depicted in Fig. (Tfa), 
"Coulombian" because of its pion propagator. The ex- 
plicit expression for the amplitude is 



iM 



where q is the momentum transfer between the nucleons, 
the Ti are usual Pauli matrices acting in the isospin space 
of nucleon i and f-^ the pion decay constant. Throughout 
this work, we employ = 92.4 MeV. The small latin 
letters refer to isospin indices of the pions as given in the 
figure. 

The amplitudes of Figs.[TJb) andjUc) are individually 
dependent on the parametrization of the pion field. The 
sum of both, however, is independent of this choice as it 
should. Therefore, we will only show results for the sum 
of both contributions for which the amplitude reads 
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(cTi • q) ((T2 • q) 



r2)}- (7) 
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Fig. 1. Leading two-nucleon contributions to n-A scattering. 
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Fig. 2. Triple scattering two-nucleon diagram. 

Here, we additionally encounter cr^, the Pauli matrices 
acting in the spin space of nucleon i. Since the propaga- 
tors contain an additional pion mass, we will refer to this 
contribution as "non-Coulombian" . 

In addition, in Ref. [20), the diagram shown in Fig. [5] 
was identified as numerically important, although in Wein- 
berg counting it appears only at order N^LO as compared 
to the leading order diagrams of Fig. [TJ This observation 
was taken as a support for the modified power counting, 
Q-counting, as discussed in the previous section. As we 
will show in section 15.31 a part of this amplitude is in- 
deed enhanced, however, not parametrically but as a con- 
sequence of the special topology of the diagram. Therefore 
we also include this triple scattering diagram, although 
this is not supported by a power counting yet. Straight- 
forward evaluation of its most dominant part gives 



iM 



-2i 



1 



(8) 



For a discussion of the full expression for the amplitude we 
refer to Sec. 15.31 It is important to note that, although the 
triple scattering diagram is enhanced compared to what is 
expected from dimensional analysis, this enhancement is 
not sufficient to fully overcome the parametric suppression 
as provided by chiral symmetry. Comparing the double 
scattering contribution, Eq. to the triple scattering 
contribution, Eq. ([8]), we still find a relative suppressiorfl 
of order 



m^/(16/^) ^ 0.14 



^ For this comparison we evaluated both operators 
isoscalar nucleon pairs and under the assumption that q ^ i 
in line with Weinberg counting. 



for 



which is close to the ratio found from the exact evalu- 
ation of the matrix elements. Thus, there is a suppres- 
sion by an order of magnitude between the two contribu- 
tions. We conjecture that the related diagram with four 
ttN interactions will be suppressed by yet another order of 
magnitude and is therefore irrelevant for this study. Note, 
however, in case of iC-nucleus scattering a resummation 
of the multiple scattering series is unavoidable [iOlHTlH^ . 
Also in Ref. [H] , it was argued that a resummation needs 
to be done including nucleon recoil effects in the propaga- 
tors, for the nucleon recoil in multiple scattering diagrams 
induces potentially important corrections of the order of 
^Jm^JniM , with being the kaon mass, as compared 
to the static contributions. 

As was argued above, all diagrams that contribute to 
lower orders than the first counter term should be included 
in this study. In addition to the diagrams discussed so far, 
there are potentially also those that are enhanced due to 
the presence of either two-nucleon or TrNN-cuts. More- 
over, there are also contributions from the Z\-isobar. The 
effect of the TrNN-cuts was discussed in detail in Ref. [43] , 
see also Ref. j44j for a related discussion. It was shown in 
Ref. [33] that, although formally enhanced, the contribu- 
tions of the TrNN-cuts to the tt-^H scattering length nearly 
vanish as a result of a cancellation of different two-nucleon 
diagrams enforced by the Pauli principle. The same argu- 
ment also applies to the scattering off ''He. In Ref. [H] 
and [1^ it was shown that diagrams with an NN-cut, of- 
ten called dispersive terms [3, and those with an interme- 
diate A start to contribute at the same order. Although 
individually sizable, the two groups of diagrams cancel 
each other nearly exactly. This cancellation should also 
appear in heavier nuclei whenever pions get scattered off 
an isoscalar NN pair. The analogous diagrams for isovector 
NN pairs were not studied yet, however, their contribution 
is intimately connected to the reaction pp — >■ ppn^. The 
total cross section near threshold in this pion production 
channel, which is the relevant quantity to estimate the 
contributions from the NN-cuts, is more than an order of 
magnitude smaller than that of pp pnir'^ [45) . In addi- 
tion, the Z\-isobar is a lot less important in that channel 
compared to the deuteron channel [35] . In total, we there- 



^ Note, the group of so-called disconnected diagrams (where 
the pion is absorbed on one nucleon and emitted from another 
one) are included in the dispersive corrections. 
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fore expect a negligible contribution from intermediate Z\'s 
and NN-cuts to pion-nucleus scattering. 

In Ref. also corrections that arise through the 
boost of the ttN interaction were identified as numerically 
significant. However, in a formalism that includes the A 
as explicit degree of freedom this is no longer the case [12] 
— see also discussion at the end of Sec. [2] 

Therefore, we may assume that the contributions listed 
in Eqs. ([5]) to are all two-nucleon contributions that 
need to be included in this calculation. 

With the help of Feynman diagrams, we have obtained 
the amplitude on two free nucleons. Following Weinberg, 
the amplitude resulting for nucleons bound can be ob- 
tained calculating the expectation value of the free ampli- 
tude using the A-nucleon bound state wave function. E.g. 
for A = 3, the wave function ipa iPi2,P3) in momentum 
space depends on the two momenta pi2 and ps, which are 
Jacobi momenta defined by 



Pl2 = 




-k2) 




PU = 




-K) 




P3 = 


3 


- ^ (fcl + fc2) 




P3 = 


3^3- 


-\{k[ + K) ■ 


(9) 



Here ki and fc/ are the momenta of the nucleons before and 
after scattering of the pion. a denotes an index labeling 
the different spin and isospin states of the three nucleons. 

For the two-nucleon contribution, making use of the 
conservation of momentum of the third nucleon, we obtain 
in this way 



-nucieonx 



(27r)3 



d^Pi2 d^P3 



X C' (P^2>P3) A^a'a (Pl'2,Pl2) V'a (Pl2,P3) (10) 

where the amplitude has also been expressed in terms of 
Jacobi momenta pi2 using the definitions of Eq. ©. We 
note that the factors {2tt)^ take into account that we nor- 
malize momentum eigenstates to S^{p — p') {{^tt)^ 5^{p — 
p')) for the wave functions (operators) throughout this 
work. Eq. ([TU)) is easily generalized to A = 4 and higher. 

The corresponding contribution to the scattering length 
can in general be determined via 



-nucLeon 




N- 



nucLeonx 



(11) 

The binomial coefficient in front is introduced to take into 
account the number of nucleon-7V-tupel in an j4-nucleon 
system. This factor is three for two-nucleon operators in 
3He. 



3.2 Leading three-nucleon contributions 

One aim of this work is to study the relative importance 
of one-, two-, three-nucleon, etc. contributions explicitly. 



91 



i <23 



Fig. 3. Coulombian three-nucleon contribution to n-A scat- 
tering. 



Therefore, we now look at three-nucleon contributions to 
pion-nucleus scattering. Based on Weinberg's original power 
counting, it is easy to identify the leading ones, which 
are shown in Figs. [3] and [M] Note that we have omit- 
ted diagrams that vanish in leading order due to thresh- 
old kinematics. Naively, these diagrams are suppressed by 
(ra^^ /rajsiY compared to the leading two-nucleon ones. A 
closer look reveals that the diagram of Fig. [3] is individu- 
ally independent of the parametrization of the pion field. 
The resulting amplitude reads 



1 



9 9 

Qi 93 



X {e'^''^[(xi.X2)T3^-K(r2-X3)rf]} (12) 

Again, the propagators resemble Coulombian ones. We 
therefore will refer to this amplitude as "Coulombian" 
three-nucleon contribution. The definition of the momen- 
tum transfers and pion isospin indices can be read off 
from the figure. Note that due to the isovector structure 
of this operator it will not contribute to pion scattering 
on isoscalar nuclei. 

Additionally to this Coulombian contribution, we found 
a set of seven further diagrams, which are individually de- 
pendent on the parametrization of the pion field. The sum 
of these diagrams is however parametrization-independent 
as is outlined in Appendix \^ and the explicit expressions 
are given in Appendix |B1 They all have in common that 
at most one of the propagators is Coulombian. In the fol- 
lowing, we will therefore refer to these diagrams as half- 
Coulombian. 

For the three-nucleon contribution, the momentum of 
the third nucleon is not conserved anymore. Therefore, 
the calculation of the expectation value with respect to 
the A = 3 wave function reads 



.3—nucleon\ 



d^p'i2 d^P3 



d Pi2 d P3 



- (2^)3 (2^)3 

X C' (Pl'2>P3)-^a'a (p{2P3>Pl2P3)'0a (Pl2,P3) 



(13) 
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4—7iucleon\ 



c 1 qi 

4 "^ 



d \ q-2 



d^p[^ Sp'^ d\'^ 



d^Pi2 d-^p3 d-^Qi 



(27r)3 (27r)3 (27r)3 

X Ipa (Pl2,P3,94) • (15) 

Here, we express the "^He wave function i/j in terms of 
Jacobi momenta. pi2 and ps are the same as in the three- 
nucleon case. 54 is the relative momentum of nucleon 4 
and the cluster of nucleons 1, 2 and 3. For details of the 
representation of the wave functions, we refer to |36) . 



Fig. 4. Coulombian four-nucleon contribution to tt-A scatter- 
ing. 



Again, the amplitude has been expressed in terms of Ja- 
cobi momenta. The contribution to the scattering length 
due to the three-nucleon diagrams is then again found by 
inserting Eq. 1^ into Eq. PT|) . 



3.3 Leading four-nucleon contribution 

In order to also confirm our conclusions on the relative 
importance of higher-body tt-A amplitudes, we also need 
a four-nucleon contribution which, of course, can only be 
relevant for scattering on nuclei with A> A. The probably 
most important contribution is Coulombian and shown 
in Fig. |4l Based on Weinberg's original power counting, 
this contribution is suppressed by [m-^/mMY compared 
to the leading two-nucleon contributions. Again, it is in- 
dependent of the parametrization of the pion field. Note 
that this diagram is not a complete set of leading four- 
nucleon terms. We will only evaluate its contribution to 
have an estimate of possible higher order contributions to 
TT-A scattering. 

It is an easy exercise to find the corresponding ampli- 
tude 



iM 



{2(5'^^ (t 
+ (-ri 



1 



? ? 2 

9i 92 93 



1 • T2) (1-3 



-{ri 



^4) (1 



2 ' 3 



3 ' 4 



'2) 



(14) 



This amplitude then enters the evaluation of the expecta- 
tion value analogously to Eqs. pil|) and (|13p for the two- 
and three-nucleon operators 



4 Numerical method 

In this section, we want to introduce briefly the numerical 
method used for the evaluation of the pertinent integrals. 
Whereas expectation values for can be easily obtained 
using standard methods of integration or using a partial 
wave decomposition, this becomes more and more tedious 
for more and more complex nuclei. One aim of this work 
was to establish a scheme to evaluate expectation values in 
momentum space based on Monte Carlo (MC) integration 
without performing a partial wave decomposition. 

In this way, we were able to generate the numerical 
expressions of the amplitude reliably with the help of a 
Mathematica script. The resulting FORTRAN code lines 
could be included in a FORTRAN code evaluating the 
high dimensional integrals given above. 

For the evaluation, it turned out that simple MC in- 
tegration requires a large number of trial points to ob- 
tain acceptably small standard deviations. In our tests, 
we found that most trial points consisted of momenta for 
which the wave functions are nearly zero. Therefore, ob- 
viously, an importance sampling similar to the Metropolis 
algorithm 47, is required to keep the computational needs 
small and increase the accuracy. 

Usually such an importance sampling is guided by the 
square of the wave function. In configuration space, this 
quantity is perfectly suited as a weight function for the 
Metropolis algorithm since the weight function is then au- 
tomatically normalized to one at least as long as the opera- 
tors are local. For momentum space, the structure is more 
complicated, since the integrals require weight functions 
with higher dimensionality as in configuration space. This 
implies that a simple square of the wave function is not 
useful for the importance sampling anymore. This prob- 
lem could be solved by performing part of the integrals 
using standard methods as has been successfully done in 
P5] . We found this approach less practical in our case, 
since the three- and four-nucleon operators would require 
to perform high dimensional integrations using standard 
integration methods. 

Our solution was to give up weight functions based on 
the wave functions of the system, but choose a rational 
ansatz instead. The parameters of the ansatz were then 
adjusted so that the standard deviation in test cases was 
minimized. In this way, we were able to improve the ac- 
curacy sufficiently. At the same time, the weight function 
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Table 1. Values for the parameters of the weight function for 
various cutoffs of the leading order chiral interaction, AVIS 
and CD-Bonn wave functions. We always chose Cp^j = 
and Cp^ = Cp' 



yl[fm-^] 


r 


C„,[fm-i] 


Cp3[fm-il 


2.0 


7.6 


1.25 


1.0 


3.0 


7.6 


1.25 


1.0 


4.0 


7.4 


2.75 


1.25 


5.0 


7.4 


2.75 


1.25 


10.0 


8.5 


5.2 


3.25 


20.0 


9.7 


S.75 


6.0 


AVIS 


7.6 


3.0 


1.25 


CD-Bonn 


7.6 


3.0 


1.25 



could be analytically normalized to one so that the calcu- 
lations became feasible. 

E.g. we choose for the importance sampling for inte- 
grals of the form of Eq. (jl3p a weight function depending 
on the four integration variables Pi — pi2 , Ps , Pi'2 a-nd p^ 

w{p[2,P3,Pl2,P3) = w{p[2,p'3,Pl2,P3) 

Stt {p^ + CpJ 

(16) 

For simplicity, the ansatz only depends on the magnitude 
of the momenta. With the parameters Cp. and r the shape 
of the weight functions can be influenced. The ansatz guar- 
antees (for large enough r) that the weight function is 
normalized to 

J dVl2C«M J rfV2d^P3w(Pl2,P3.Pl2,P3) = 1- (17) 

A detailed description of our tests of this method can 
be found in |49j . Here we only summarize the resulting pa- 
rameters in Tabled] For the more simple integral Eq. ([T0|) . 
the weight function was simplified in the obvious way just 
dropping terms depending on Pg. We used the same pa- 
rameters in both cases. 

Still, the size of the integrand is driven by the size of 
the wave functions. This is reflected in a strong depen- 
dence of the parameters on the wave function used. In the 
table, we give our choice for the leading order chiral wave 
functions with different cutoffs A and the wave functions 
obtained based on AV18 and CD-Bonn. For the higher or- 
der chiral wave functions, we used the same parameters 
as for the leading order A = 2 fm~^ wave functions, since 
the momentum dependence is similar in this case. 

As a first test, we compare results for obtained with 
a traditional partial wave (PW) decomposition and with 
the MC method. The partial wave decomposed amplitudes 
are listed in Appendix [C] In Table [21 the results are com- 
pared. One can see that the MC results agree well with 
the PW values and that these agree well with the previ- 
ously obtained values of ^0]. We note that 2/3 of the MC 
results are within Icr of the PW result. We found a few 
values with more than 2a deviation. These outliers are 



generally expected for a MC calculation and indicate that 
the probability distribution is not normal, but has more 
extended shoulders. Note, as usual the numerical accuracy 
of the MC calculation can easily be enhanced by increas- 
ing the number of runs — we stopped our evaluations as 
soon as the numerical accuracy was higher than the the- 
oretical accuracy of the calculation, which is of the order 
of a few percent as discussed below. 



5 Testing the counting schemes 

5.1 Cutoff dependence and estimate of the leading 
counter term 

As was argued in the introduction, from the regulator de- 
pendence of the contributions from the various diagrams 
it is possible to estimate the size of the leading 4N7r7r 
counter term contribution. Since the estimated contribu- 
tion of this term differs drastically between Q-counting 
and Weinberg counting, a study of the cutoff dependence 
provides a non-trivial test of the counting schemes. 
The leading chiral NN interaction is given by 

In order to obtain a meaningful Schrodinger equation, it 
is necessary to introduce a regulator. We here perform 
regularization by a smooth cutoff function given by 

/(p)=exp(^-g)') (19) 

which depends on a cutoff parameter A. The potential is 
then replaced by 

V{p',p)^ f{p')V{p',p) f{p) (20) 

In this section, we used gA — 1.2834, ttItt = 139.57 MeV, 
and /tt = 92.4 MeV for the leading order potential and 
the pion scattering amplitudes which is close to the pa- 
rameters used in [5D]. Note that we use slightly different 
parameters for studying '^He and ^He below. 

Since we are only interested in the ^Si-'^Di partial 
wave, the deuteron channel, we arbitrarily set Ct ~ 
and fit Cs so that the deuteron binding energy is fixed 
to a experimental value for a given value of the cutoff 
parameter A. The fit results are given in Table [31 

We are now in the position to study the A dependence 
of the scattering length contributions due to Eqs. (O, (0), 
and ([5]) based on the LO potential. For this, we have ad- 
justed the potential so that the deuteron binding energy 
is close to the experimental one. For the Coulombian di- 
agram, the cutoff dependence has already been studied 
in [50| in momentum space. This diagram and the triple 
scattering diagram, as given in Fig. [51 have also been in- 
vestigated in Refs. [STIISS] with the result that, for wave 
functions based on the l7r-exchange interaction, the re- 
sults become independent of the cutoff in the limit of large 



S. Liebig et al.: Towards a high precision calculation for the pion-nucleus scattering lengths 



9 



Table 2. Comparison of PW and MC results. The scattering length contribution of the different two-nucleon operators Figs[l] 
(a), (b) + (c) andOis compared for different cutoffs and phenomenological wave functions. 



A[im-^] 


m 


[10-3 m-i] 




[10-3 m-i] 


J?7V 

12} 


[10-3 m-i] 




PW 


MC 


PW 


MC 


PW 


MC 


3.0 


-21.37 


-21.16(22) 


-0.666 


-0.6658(4) 


2.43 


2.433(3) 


4.0 


-20.02 


-19.65(9) 


-0.902 


-0.9017(7) 


1.77 


1.767(2) 


5.0 


-19.75 


-20.13(35) 


-0.889 


-0.8904(6) 


1.61 


1.613(3) 


10.0 


-20.72 


-22.68(136) 


-0.824 


-0.8232(14) 


2.35 


2.345(7) 


20.0 


-20.78 


-20.60(52) 


-0.867 


-0.8729(40) 


2.49 


2.459(29) 


AV18 


-19.62 


-19.45(13) 


-0.749 


-0.7493(6) 


1.63 


1.631(3) 


Nijm93 


-19.84 


-19.40(11) 


-0.743 


-0.7425(3) 


1.72 


1.724(3) 


CD-Bonn 


-20.20 


-19.99(8) 


-0.553 


-0.5521(4) 


1.92 


1.924(2) 


CD-Bonn [20^ 


-20.20 




-0.55 









-19.6 




2 4 6 8 10 12 14 16 18 20 

A [fm-i] 



Fig. 5. Contribution 1^23 "^^^ ^'^ amplitude Eq. © to the 
TT-^H scattering length depending on the cutoff A of the LO 
potential. The result for the Nijmegen 93 is also shown "(+)"• 
ft is independent of A and is positioned arbitrarily on the plot. 



2.6 




A [fm-^] 

Fig. 7. Contribution a^2H '^^^ 1° amplitude Eq. ^ to the 
TT-^H scattering length depending on the cutoff A of the LO 
potential. The result for the Nijmegen 93 is also shown "(-I-)". 



-0.65 




Fig. 6. Contribution fflJl2H "^^^ amplitude Eq. ([7| to the 
TT-^H scattering length depending on the cutoff A of the LO 
potential. The result for the Nijmegen 93 is also shown "(+)"• 

A. In Figs.jSJini and [71 this result is seen once more. Note 
that the scales are entirely different in these figures. Ad- 
ditionally, we observe that the same is true for the con- 
tribution due to Eq. ([T]). We observe a mild oscillation of 
the result, the amplitude of which is approximately 5 % 
of the respective contribution — 1 x IQ-"^ in absolute 



terms. From the comparison in Table [51 we observe that 
the leading order results are in good agreement with the 
previous potential model results. This is also true for the 
amplitude in Eq. From this we conclude that a theo- 
retical accuracy of 1 x 10^^ can be reached at most 
from a study of pion-nucleus scattering, which means that 
a 5% accuracy can be reached. 

As was already argued in Sec.[2l the dependence of the 
few-body operators on the regulator used for the deuteron 
wave function can give us an idea on the numerical size of 
the leading counter term contribution. In Weinberg count- 
ing the leading isoscalar counter term is expected to be 
suppressed by a factor (rriT^/mjsiY ~ 2% compared to the 
numerically leading double scattering term. This is fully 
in line with the observed amount of A dependence, as de- 
scribed in the previous paragraph, but in gross disagree- 
ment to the expectations of Q-counting, where a counter 
term contribution of the order of 0.2 % would be expected. 
Thus, from the study of the cutoff dependence we con- 
clude that Weinberg counting provides the more accurate 
estimate for the leading counter term contribution. 
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Table 3. Values for Cs depending on the cutoff A. For these 
fits, the experimental binding energy Ed ~ 2.225 MeV is re- 
produced. We fixed Ct = in all cases. 



A[im-^] 




3.0 


-34.2225 


4.0 


48.1751 


5.0 


562.089 


10.0 


-50.9683 


20.0 


-92.1179 



5.2 Dependence of the tt-^H scattering length on the 
deuteron binding energy 

Following up on the discussion in Sec.[5J we now study the 
energy dependence of the ratios of few-nucleon contribu- 
tions. If Q-counting were operative, the relative scaling of 
Diagrams Figs.[TJa) and[ljb)4nic) would strongly depend 
on the binding momentum. A straightforward analysis of 
Eqs. © and ([7]) shows that Diagram (a) of Fig. [1] scales 
as 

(la) 2 ; 2 

while the sum of Diagram (6) and (c) scales as 

(Ibc) 2 / 2 

Thus, for the binding energy dependence of the ratio of 
these two contributions we find 

a^^A^ oc / (^^^"^^^"^g) (21) 

ai';^ ^ I (Q-counting) ' 

where we used the Q-counting relation (7^7 together 
with Eq. (0) . Analogously one gets 

/ const. (Weinberg) 
<^i'A ^ (g-counting) ' ^ ' 

using the explicit expression for the contributions of the 
individual diagrams. 

We have studied the binding energy dependence of 
both classes of diagrams based on the leading chiral NN in- 
teraction. To this aim, we have adjusted the contact inter- 
action acting in the deuteron channel so that the deuteron 
was bound with unphysically large and small binding en- 
ergy (see Table E]) . It was then an easy exercise to cal- 
culate the contributions of Figs. [H^a), [IIb)-|-[T](c) and [2] 
to the scattering length depending on the binding energy 
and explicitly compare to the expectations from Q~ and 
Weinberg counting. 

Based on this observation, we are now in the position 
to look in more detail at the dependence of the various 
contributions on the binding energy. For this, we arbitrar- 
ily choose A = 20 fm~^, which is in the region where the 
results are almost independent of the cutoff. We start dis- 
playing the binding energy dependence of the individual 
contributions to the scattering length in Fig. HI We ob- 
serve that, independent of the binding energy, the contri- 
bution of the Coulombian diagram is the most important 



Table 4. Values for Cs depending on the chosen binding en- 
ergy of the deuteron Ed. For these fits, we fixed Ct ~ and 
yl = 20 fm-^ 



ii£j[Me VJ 


Os[<oceV ] 


0.002 


-71.7689 


0.005 


-72.1400 


0.01 


-72.5574 


0.02 


-73.1460 


0.05 


-74.3088 


0.1 


-75.6126 


0.2 


-77.4471 


0.5 


-81.0703 


1.0 


-85.1564 


5.0 


-103.225 


10.0 


-118.793 


20.0 


-147.032 


30.0 


-178.149 


40.0 


-217.147 


50.0 


-271.056 



two-nucleon contribution. In contrast to the naive power 
counting estimates, the amplitude of Eq. (|5]) — the triple 
scattering diagram, depicted in Fig. [5] — is the next im- 
portant one. Still, it is suppressed by one order of mag- 
nitude compared to the leading Coulombian two-nucleon 
diagram. The contribution of this diagram will be dis- 
cussed in detail in Sec. 15.31 The amplitude Eq. (O — 
shown in Fig. [ljb)-|-[llc) — gives an extraordinarily small 
shift of the scattering length. From the observation that 
this suppression is not strongly depending on the bind- 
ing energy, we conclude that this suppression is unrelated 
to the binding momentum as suggested in Ref. [3D], but 
probably accidental. 

In order to be more quantitative on the relative sup- 
pression of these contributions, we show in Fig. [HKd) the 
ratios of the shifts of the scattering lengths due to Eqs. (I?]) 
and ([5]) and the one due to Eq. ^ . Based on Q-counting, 
the ratio for Eq. (O should scale like Ej (c.f. Eqs. for 
small energies, the one for Eq. ([5]) should scale like y/Ed 
(c.f. Eqs. (l22t ). Our explicit calculation, however, shows 
a much weaker dependence on the binding energy. This 
is in strong contrast to the expectation from Q-counting. 
Therefore, we conclude that Q-counting is not realized for 
pion scattering on light nuclei, i.e. ^H. 

It is instructive to investigate further the source of 
the residual binding energy dependence reported above. 
From this we will find that the ratios of contributions to 
the scattering length has a logarithmic and not a power 
law dependence. In addition, we will be able to show that 
the physical deuteron binding energy is already beyond 
the range of applicability of heavy pion effective field the- 
ory. To this aim, we study the ratios within an analyti- 
cal model of the deuteron wave function. To get there we 
start from pionless EFT to obtain the deuteron based on a 
zero range approximation of the NN interaction. However, 
from different investigations it became clear that the bind- 
ing momentum cannot be the only scale affecting the wave 
functions of the deuteron in a model independent way: tt- 
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Fig. 8. Contributions a^l-^ (a), af^s^ (b), and cl^2ii (c) due to the amplitudes of Eqs. ((B]), 0, and ([SJ to the tt-^H 
scattering length depending on the binding energy Ed of ^H. The results for Nijmegen 93 are also shown "(+)"• (d) shows the 
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Fig. 9. Comparison of the Hulthen deuteron wave function to the S-wave part of several LO wave functions at the physical 
deuteron binding energy (a) and at Ed = 0.01 MeV (b). The Hulthen wave function is shown for /? = 1.7m-n. We also show the 
deuteron wave functions based on contact interactions (/? — >■ oo). 



scattering at threshold was studied based on the per- 
turbative treatment of pions in [53] and using heavy pion 
EFT in Refs. [5¥lf55] with the result that a counter term is 
required in leading order where two-nucleon diagrams con- 
tribute. In practice, this would imply that an extraction 
of a*^+^ is unfeasible based on tt-^H atoms. However, it was 
realized that this problem is tamed once pions are treated 
non-perturbatively |50II51[I521 I56'. In this case, the short 
distance wave function is affected by l7r-exchange in 
such a way, that counter terms are not required to obtain 
cutoff independent results. No contradiction to the naive 



power counting of Weinberg is observed. Obviously, the 
pion introduced scales into the wave function beyond the 
binding momentum. Therefore, we add a range factor to 
the vertex function, so that we are able to introduce the 
intrinsic non-perturbative scale that enters through pion 
exchange. In momentum space, the wave function then 
reads 



np)=N{j,l3) , (23) 
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where the normalization factor N is fixed to 

7V(7,/3)' = 8^7/3(7 + /?)' 

by the normahzation condition for the deuteron wave func- 
tion. This is a wave function of the Hulthen type ^5T . In 
Refs. [56.,50.,58j . it was shown that for radii larger than 
0.6 fm the LO chiral wave functions are basically inde- 
pendent of the regulator used for their construction. We 
therefore fix /? by fitting to the tail of the wave function 
at Ed = 0.01 MeV. This way we find 



/3 = 1.7m7r 



(24) 



It is reassuring that /3 turns out to be of the order 
of the pion mass. If the proposed picture is correct, the 
scale /3 should be independent of the deuteron binding en- 
ergy. We confirmed that this is indeed the case as long 
as we do not go to very large binding energies above the 
physical one. Fig.[H]shows the S'-wave deuteron wave func- 
tions compared to the Hulthen ansatz. For both binding 
energies, we use the same /3. As one can see, the Hulthen 
and chiral wave functions nicely agree for larger distances. 
In Fig. EJa) one can also see that this long range part is 
independent of the cutoff used. 

The limit of a point like vertex is achieved by the limit 
/3 — oo. In this limit the wave function of Eq. (1^51) be- 
comes identical to the one used in theories with perturba- 
tive pions [53] as well as the one used when treating the 
pion as heavy field [M1I55) . For completeness, we also show 
this wave function in Fig. M,^) ■ It is obvious that such a 
simplistic wave function is not a good approximation to 
complete chiral wave functions. The most important ef- 
fect is a reduction of the long range part, which becomes 
necessary to insure the correct normalization of the wave 
function. 

The wave function of Eq. (l23l) is still sufficiently simple 
that analytic calculations can be performed for the various 
matrix elements discussed in this section. In particular we 
find 



(i'^) (1 + ^) 



In 



4x 



(25) 



where we introduced the dimensionless parameter x 
j/P and 



K = 



1 



87r2 (1 + m^/2mN) 
In addition, we find 

(1 + ^) 9a 



(Ibc) _ 
7r — 



X On 



(l-a;)2 12 
/ {1 + X + ffi^f 
\ (2x + m7r)(2 + TOjr) 

m7r(l — a;)2 
{1 + X + TO7r)(2a; -I- m7r)(2 -I- m^) 

with ttItt = mj^/P and 

„(2) (1+^) 



,(26) 



(1 



47r/2 



In 



1 



:ln 



2x 
x + l 



(27) 



We used these analytic results, found with the Hulthen 
wave functions, to predict the ratio 



-2H 



(28) 



In Fig. [Tol a'), this analytical result is compared to a LO 
calculation for which we neglect the D-wave contribution 
of the wave function. Although the results do not agree 
quantitatively, it is obvious that the simplified calcula- 
tion based on the Hulthen ansatz is able to describe the 
energy dependence qualitatively. The inclusion of the D- 
wave changes the result for Diagram[TJa) only marginally, 
however, the full contribution for Diagram[TIb)-f{TIc) even 
changes its sign. Let us now focus on the contributions 
from only the deuteron S'-wave. From Eqs. ([25]) and it 
follows directly that the suppression of Diagram[TJb)+mc) 
compared to the Coulombian [Ija) is only logarithmic, in 
line with Weinberg counting, and not power law (7^) as 
predicted by Q-counting. Close inspection reveals that the 
bulk of the suppression of ((Tt»c) with respect to ((TJi) comes 
from spin-isospin factors leading to the factor of 1/12 

Such kind of accidental suppression, a power 



Ibc) 



m a 

counting can not capture. 

It is also very interesting to investigate various limits of 
Eqs. and As mentioned above the expressions for 
the wave functions relevant for theories with perturbative 
pions, as used e.g. in Ref. [53], are recovered when taking 
the limit /3 — 00. Then the above expressions collaps to 



K, X ln(4a;) , 



(29) 



and 



(Ibc) 



.9\ 



12 

X I In (2(20: -hm^)) + 



{2x + m^) 



(30) 



Identifying /3/2 with A*, the given equations agree with 
those of Ref. [53] . In addition we find 



,(2) 



^H,LO 



47r/2 



In (2) 



Contrary to the theory with perturbative pions, in the 
heavy pion effective field theory the Diagrams [ljb)-|-[l](c) 
do not appear explicitly but are absorbed in local counter 
terms. A comparison of the given leading order expressions 
with the full result reveals that, for Diagrams [TJb)-fnic), 
there is no regime of binding energies where the given for- 
mula well represents the full result. Contrary to this, the 
leading order expression for Diagram [TJa) works well for 
very small values of the binding energy (see Fig. [TUF b)). 
This is to be expected since there must be a kinematical 
regime where the heavy pion effective field theory is appli- 
cable. However, we find that the expressions break down 
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Fig. 10. Comparison of the ratio r 
and the S-wave part of the LO wave function for A 



(a) and of t[^2H C^) obtained from the Hulthen deuteron wave function 
20 fm~^ depending on the deuteron binding energy Ed- The Hulthen 



result is shown for /3 — 1.7m,r- In (b) the results are also compared to the expansion of Eq. (|29p . 



already for very low values of the binding energy. Already 
for a binding energy of 0.1 MeV the leading order expres- 
sions show a significant deviation from the full result. We 
therefore conclude that for accurate calculations of tt-^H 
scattering neither the treatment with perturbative pions 
nor heavy pion effective field theory are applicable for the 
physical deuteron. 



5.3 Analysis of the triple scattering diagram 

The observation that the diagram shown in Fig. [5] is sig- 
nificantly enhanced compared to the expectation based on 
Weinberg counting was taken as a further support for Q~ 
counting [50] . In this subsection, we will demonstrate that 
the triple scattering diagram is not enhanced as a result 
of the smallness of the deuteron binding momentum but 
because of the special topology of its loop diagram. 

The loop that appears in the triple scattering diagram 
may be written as 



1 



(v ■ I ~ uj — ie) 



1 



(m^ — P — ie)(m^ ~ (J ~ QY ~ 



(31) 



where v — (1,0,0,0). In Ref. [5S] a general solution for 
this integral is given. In the kinematics relevant for -n-A 
scattering at threshold we find 



with 
6h 



/o(m^, 0, -q^) = —— + SIq, 



— — — r ( dx I arctan | — — : 

8^2|q| y \smix)\q\ 



(32) 



(33) 



In Ref. '20.. only the first term on the right hand side was 
included. Note that in the limit of heavy pions, as used 
in Ref. [5S], the contribution of SIq vanishes. Dimensional 
analysis, which is the basis of Weinberg counting, allows 



us to estimate integrals. In case of Iq this analysis gives 
(assuming q ^ niTr) 



(34) 



Jo - l/((4^)'m,) 



Clearly, the first term on the right hand side of Eq. 
is enhanced by a factor of 27r^ compared to this estimate. 
The remainder, SIq, on the other hand is numerically fully 
in line with the estimate p4p . The power counting can 
only capture parametric suppressions. We therefore con- 
clude that the fact that SIq behaves in accordance with 
Weinbergs power counting is a further support for its ap- 
plicability. However, in a high accuracy calculation for 
pion-nucleus scattering the first term of Eq. (j32p is to 
be kept and it is this piece that we have in the list of 
operators to be included in the calculation — c.f. Eq. ([8]). 

How can we understand the large enhancement of a 
part of Jo? To see this, we first observe that the enhanced 
part of Jo can be directly calculated from the full Feynman 
integral by only keeping the term that corresponds to the 
nucleon pole. For this piece we get 



xnucl. pole 
^0 




(35) 



The implicit assumption behind any dimensional analysis 
is that integrals, once converted into dimensionless expres- 
sions, are of order 1. And indeed, if the given integral were 
of order 1, the full expression would be perfectly in line 
with power counting. However, one finds 




We trace the appearance of this large result to the presence 
of an integrable singularity at a: = 1 in the expression 
given above. Indeed, 80 % of the exact result originate 
from values of x < 2. 

It is interesting to note that enhancements by fac- 
tors of TT were already observed to emerge also in pion 
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Table 5. Values for Cs and Ct depending on the cutoff A of 
the LO potential. 



yl[fm-i] 


Cs[GeV-2] 


CT[GeV-2] 


2.0 


-83.6941 


2.63787 


3.0 


-29.0931 


16.3942 


4.0 


86.3303 


52.6427 


5.0 


-435.354 


-122.611 


10.0 


-39.8356 


6.36715 


20.0 


-66.2861 


-4.03158 



loop contributions to the NN potential [60] , the scalar nu- 
cleon form-factor [61j and the tt" photoproduction ampli- 
tude from similar topologies as those discussed here. 
A deeper understanding, when these dimensionless factors 
appear, would be very desireable. 

It is also important to note that integrals of the same 
topology are also relevant for the reaction NN — >■ NNtt, 
though in different kinematics. For a discussion of how to 
treat these large momentum transfer reactions in ChPT, 
see Refs. [551H5] . Loops for this reaction were studied in 
Refs. [W, 64,65,6^. It is conceivable that the problems in 
understanding quantitatively especially the reaction pp — >■ 
ppiT^ are connected to the same enhancement of loops as 
discussed in this section. 



6 Few-nucleon wave functions 

Before we can give results for the scattering lengths, we 
need to specify the input going into our calculations of the 
few-nucleon wave functions in more detail. We will again 
study the results for LO wave functions. In contrast to ^H, 
we restrict to interactions in "^Si-^Di and ^Sq partial waves 
requiring fits of the strength of both contact interactions, 
Cs and Ct- We determined these constants fitting to the 
binding energy of and the ^Sq neutron-proton phase 
shift at a laboratory energy of 1 MeV. The potential was 
given already in Eq. (|18p . For the new fits, we used for 
historical reasons a different, more sharp regulator 

/(P)=exp(^-g)') , (36) 

gA — 1.29 (utilizing the Goldberger-Treiman relation jB71) 
and = 138.0 MeV. The results for Cs and Ct are 
summarized in Table [5j 

Based on these fits, it is a straightforward task to cal- 
culate the binding energies of "^He and ^He. To this aim, 
we have solved Faddeev/Yakubovsky equations in momen- 
tum space following [36 . Thereby, we subtracted the poles 
of unphysical spurious NN bound states from the 2N t- 
matrix as outlined in |37| . As was already shown in the 
same reference, the binding energies become rather inde- 
pendent of A for large A. However, in LO, the binding 
energies of '^He and ^He do not well reproduce the ex- 
perimental values. Table |6] shows our results for the var- 
ious NN potentials used. In LO, the '^He binding energy 
is varying for the different cutoffs by almost 6 MeV. Such 



Table 6. Summary of the ^He and *He binding energy results 
for the LO, NLO, and N^LO chiral interactions, AV18 and CD- 
Bonn. For the LO interaction, the cutoff A is given in [fm~^]. 
For the chiral interaction, the Lippmann-Schwinger cutoff A 
and spectral function cutoff A is given in MeV The binding 
energies are given m MeV. For ^He, we have not performed 
calculations for all cutoffs in LO. 





A / A 


B(=^He) 


BC'He) 


LO 


2.0 / - 


11.042 


39.88 


LO 


3.0 / - 


6.878 


20.25 


LO 


4.0 / - 


6.068 


17.08 


LO 


5.0 / - 


5.987 


16.48 


LO 


10.0 / - 


5.611 


15.05 


LO 


20.0 / - 


5.429 




NLO 


400/500 


7.678 


28.57 


NLO 


550/500 


6.991 


24.38 


NLO 


550/600 


7.051 


24.72 


NLO 


400/700 


7.699 


28.77 


NLO 


550/700 


7.090 


24.94 


N^LO 


450/500 


7.717 


28.04 


N^LO 


600/500 


7.740 


28.11 


N^LO 


550/600 


7.722 


28.28 


N^LO 


450/700 


7.726 


27.65 


N^LO 


600/700 


7.808 


28.57 


CD-Bonn 




7.719 


28.28 


AV18 




7.736 


28.36 


Expt. 




7.718 


28.30 



a large variation can be expected in low orders, since the 
binding energies are specifically sensitive to changes of the 
potential 68 . Note that we did not include the Coulomb 
interaction in these LO calculations, whereas we did in- 
clude the Coulomb interaction for the other orders and 
the phenomenological NN forces. In NLO, for a smaller 
range of cutoffs, the variation is reduced, but still visi- 
ble. In NLO, again, there is an appreciable deviation from 
the experimental values. It is only in N^LO, that three- 
nucleon forces (3NF's) contribute [^155] . which ensure by 
construction that the binding energies are close to the ex- 
perimental results. For the phenomenological interactions, 
it is by now standard to augment the Hamiltonians by 
phenomenological 3NF's [70ll36l[7T] mostly based on the 
Urbana [75] or Tucson-Melbourne |73] models. Also here, 
the 3NF's are then adjusted so that the binding energies 
of^Heand ''He are close to the experimental values. Based 
on these adjustments, we are now in the position to calcu- 
late shifts of the pion-nucleus scattering length due to the 
few-nucleon contributions based on LO, NLO, and N^LO 
chiral wave functions and on state-of-the-art phenomeno- 
logical ones. 

7 Two- and three-nucleon contributions to 
TT-^He scattering 

We now study the few-nucleon contributions to tt-'^Hc 
scattering in more detail with the goal to get a better. 
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Table 7. Summary of the shifts of the Tr-'^He scattering length due to the few-nucleon corrections ^7^— ^He' '^tt— ^Hc' 

'^^aHc' '^S^Hc' ^^'-^ '^S^Hc' cutoff A is given in fm^^, for NLO and N^LO, both cutoffs (A/A) are given in MeV. 

Central values and standard deviation for the scattering length results are given in units of flO~^ ^1 • 





A/A 


&) 


etc) 

",r-3He 


l2js) 


f2lv) 


J3} 




CD Bonn 




-25.08(6) 


-0.329(1) 


2.769(2) 


0.890(1) 


-4.020(73) 


-0.789(6) 


AV18 


_ 


-24.13(9) 


-0.884(1) 


2.286(3) 


0.788(1) 


-3.536(36) 


-0.728(4) 


T 


on/ 




1 71 1 n ^ 




i.yooi^i ) 


— D. / OOi^lZo ) 


— z.ooy i^D ) 


LO 


3.0/ - 






2.445(1) 


1.553(1) 




— 1.762(4) 


LO 


4.0/ - 


-20.70(25) 


-0.327(1) 


1.077(3) 


1.581(1) 


-3.620(25) 


-1.556(15) 


LO 


5.0/ - 


-20.33(13) 


-0.120(1) 


0.868(3) 


1.812(1) 


-3.830(93) 


-1.740(15) 


LO 


10.0/ - 


-21.92(99) 


0.128(3) 


0.770(11) 


2.576(5) 


-4.283(302) 


-1.872(57) 


LO 


20.0/ - 


-19.22(146) 


0.296(18) 


-0.392(98) 


4.005(76) 


-3.835(894) 


-1.743(158) 


NLO 


400/500 


-25.36(4) 


0.828(1) 


3.117(1) 


0.984(1) 


-3.934(29) 


-0.695(1) 


NLO 


550/500 


-24.33(4) 


-0.061(1) 


2.714(1) 


0.645(1) 


-3.431(141) 


-0.374(1) 


NLO 


550/600 


-24.05(2) 


-0.037(1) 


2.637(1) 


0.671(1) 


-3.245(22) 


-0.397(4) 


NLO 


400/700 


-25.23(3) 


0.847(1) 


3.085(2) 


1.002(1) 


-3.898(21) 


-0.720(2) 


NLO 


550/700 


-24.05(5) 


-0.020(1) 


2.564(1) 


0.692(1) 


-3.311(15) 


-0.435(4) 


N^LO 


450/500 


-25.76(4) 


0.642(1) 


3.189(2) 


0.987(1) 


-3.979(7) 


-0.721(2) 


N^LO 


600/500 


-25.60(5) 


-0.021(1) 


3.039(3) 


0.779(1) 


-3.826(58) 


-0.496(4) 


N^LO 


550/600 


-25.55(2) 


0.233(1) 


3.104(3) 


0.952(1) 


-4.057(75) 


-0.708(4) 


N^LO 


450/700 


-25.25(4) 


0.611(1) 


3.104(1) 


1.052(1) 


-4.038(15) 


-0.806(1) 


N^LO 


600/700 


-25.51(7) 


0.094(1) 


3.022(2) 


0.985(1) 


-4.056(53) 


-0.734(2) 



quantitative understanding of the relative importance of 
N- and (N + l)-nucleon operators. 

The results for •^He are summarized in Table [71 To ob- 
tain these values, we have used the Monte Carlo scheme 
introduced in Sec. H) The table gives the averaged results 
together with the estimate of the standard deviation. In 
all cases, we have performed several independent runs and 
checked that the spread of the different results is in rea- 
sonable agreement with the expectations from our esti- 
mates of the standard deviation. The table distinguishes 
the shifts of the tt-'^Hc scattering due to Eqs. ^ ('^U^shg)' 
O (a|^3Ho)' ® (isoscalar part a|^^ifj(, and isovector part 
o^^He' respectively), and ([T^ (o^ajj^), and due to the 
sum of the contributions listed in Appendix [B] (a j^,^jj^). 

Based on reasonable assumptions for a^^^ the one- 
nucleon contribution to the scattering length was found to 

be a^^^3Hc = (92 ± 15) x lO-^m"! [T!]. The uncertainty 
in this result is mainly due to the uncertainty in a^"*"^ mul- 
tiplied by 3 as follows from Eq. ([TJ. On the other hand, 
it was realized in Refs.[H rHl[7Fj that the inclusion of the 
leading IV effects in tt-N scattering leads to the replace- 
ment of a^^^ by a^^^ in Eq.(IT]). In addition to this, at the 
same order there is also an isospin violating electromag- 
netic correction to tt-'^Hc scattering: — q;/2/2/(1 + g™^'" ) 
with a = 1/137 and LEG /z = (-0.97 ± 0.38)GeV-i 
[7]. The latter, however, gives a relatively small shift of 
the scattering length by (0.5 ± 0.2) x lO^^m^^. A recent 
systematic analysis of isospin violating effects in tt-N scat- 



tering up to NLO OE] resulted in updated values for a'^+' 
and a(-) [75] fro m a combined analysis of pionic hydrogen 
and deuterium data: 

= (l±l)xlO"^m^\ a^-^ = (86.5±1.2)xlO"^TO^^ . 

In the same works numerically significant subleading IV 
corrections were identified. Those may be included here 
by changing 

a(+) + and a^^^ ^ a!--^ + Ao!^-'^ , 

with zia(+) (-3.35 ±0.28) and Z\a(-) = (1.39 ± 1.33), 
both values again in units of 10^"^ m^^. Equipped with 
these numbers we get an updated value for the one-nucleon 
contribution a-^^,^^ = (88 ± 4) x lO^^ m^i. 

Also in Ref. [71], the contribution of the two-nucleon 
diagrams was estimated based on approximated wave func- 
tions for CD-Bonn. Their result is cP'^],^ = — 26x 10^^ m^^ 
which is in acceptable agreement with our full calculation. 
In this work, we aim at the theoretical improvement of the 
result of Ref. [TJ] in several aspects. First, using chiral 
nuclear wave functions up to N^LO allows us to analyze 
systematically the model dependence of our results. Sec- 
ond, the empirical enhancement of the triple scattering 
diagram, discussed in subsection 15.31 calls for an inclu- 
sion of this two-nucleon operator also in tt-'^Hc scattering. 
We find that the isoscalar part of the triple scattering 
diagram reduces the leading double scattering contribu- 
tion by about 12% which is fully in line with the corre- 
sponding contribution to tt-^H scattering. Moreover, the 
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Leading few-nucleon contributions u 

^ 7r — ^Hc 

3jj^ to the TT-^He scattering length. The crosses denote the 
results obtained by using the LO chiral wave function (triangle 
up: AV18, triangle down: CD-Bonn). 



triple scattering diagram is even further enhanced to a 
quite sizable 15% contribution once the isovector part is 
included. In addition, for the first time, we investigate the 
role of the leading three-nucleon contributions. From Ta- 
ble [7] it becomes clear that three-nucleon contributions are 
suppressed compared to the two-nucleon ones. Again, the 
Coulombian contributions ai^^ljj^,, o^^^shc ^^'^ ai^^'au^ are 



-3He 



a lot more important than the non- or half- Coulombian 
ones although the binding momentum is larger for "^He 
than for ^H. 

It was already discussed in Sec. [5] that, in case of isovec- 
tor nuclei, counter terms start to contribute from lower 
orders than in isoscalar nuclei. Due to this fact the theo- 
retical accuracy of the extraction of the ttN low energy pa- 
rameters from TT-^He is, unfortunately, significantly lower 
than from ir-^He or 7r-''He. The contribution of the isovec- 
tor counter term in 7r-^He scattering can be estimated 
using dimensional analysis to be ^ mT^/m^ 



-3Ho ~ 

4 X 10~^'m~^. On the other hand, the isovector contact 
term is expected to be enhanced by itin /rriT^ compared to 
its isoscalar counter part, which was estimated above (see 
sec. l5.ip to be of order 1 x 10~^m~^. This would provide 
us with an estimate of 7 x 10"'^ "^^^ ■ Both numbers are 
consistent and we use the latter uncertainty below. 

To support our uncertainty estimates it is important 
to study the cutoff dependence of the various scattering 
length shifts and to compare it to the counter term con- 
tribution. The cutoff dependence is shown in Figs. [TT] and 
fT2l We noticed that the results depend strongest on the 
cutoff of the Lippmann-Schwinger equation A. Therefore, 
we plotted the results depending on this cutoff for NLO 
and N^LO neglecting the mild dependence on the spectral 
function cutoff A (see [28] for more details on the definition 
of this cutoff). We also included the results for CD-Bonn 
and AV18 in these figures, which we arbitrarily positioned 
left of the other data. 

Fig. [TT] shows the most important two- and three- 
nucleon contributions, the Coulombian ones. It becomes 
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Fie. 12. Few-nucleon contributions a W., and a' „ as 
a function of A. The crosses, the triangle up and the triangle 
down mark again the LO chiral wave functions, AV18 and CD- 
Bonn, respectively. 
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Fig. 13. Results for the leading two- and three-nucleon con- 
tributions and for the isoscalar part of Fig. [2] respectively. 
Crosses: LO chiral wave functions, circles: NLO, squares: 
N^LO, triangle down: CD-Bonn, triangle up: AV18. For the 
NLO (N^LO) four (two) results were selected since they are 
very close and therefore not distinguishable. 



clear that most of the dependence on the cutoff is for cut- 
offs below 5 fm^^. The suppressed contributions for the 
non- or half- Coulombian diagrams are shown in Fig. 1121 
Also here, we observe that the cutoff dependence becomes 
mild for larger A. But the results based on the LO wave 
functions are not always in good agreement with the ones 
for phenomenological and higher order wave functions. 
Apparently, the cutoff dependence is not as strong as other 
higher order contributions for these diagrams. Obviously, 
the small size of the contributions amplifies small effects. 
Especially, the counter term contribution is significantly 
larger than the few-nucleon terms presented in FigfT^ 

Quantitatively more significant is the cutoff depen- 
dence of the Coulombian diagrams for small A. We found 
that for these contributions the cutoff dependence is driven 
by the binding energy of "^He. This is shown explicitly in 
Fig. [13] where the results of Fig. [TT] are plotted again, 
this time depending on the binding energy. Basically, the 
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smaller is the cutoff the larger is the binding energy (see 
Table [6]) and the larger is the scattering length for the 
dominant two- and three-nucleon contributions. Whereas 
this fact does not change our conclusion on the power 
counting below, it will be interesting for the extraction 
of TT-N scattering lengths from data on light nuclei. Ob- 
viously, the uncertainty due to higher order contributions 
to the wave function can be reduced by the requirement 
that the binding energy is correctly described. In this way, 
the effective dependence of the scattering length shifts on 
the order of the interaction is reduced to approximately 
1 xlO~^m~^. This uncertainty is smaller than the esti- 
mated contribution of the contact term and therefore ir- 
relevant for the estimation of total theoretical accuracy. 

The net contribution of two- and three-nucleon terms 
to TT-^He scattering length can be read off from Table [7] 
(the numbers are in 10~^m~^) 

a^f^Ho^^ = (-25.6 + 4.0 + 0.3) + (-4.0 - 0.7) ± 7 

= -26.0±7. (38) 

Here the numbers in the first (second) bracket correspond 
to the average results for two(three)-nucleon contributions 
calculated with N^LO wave functions. The uncertainty 
due to the use of different wave functions is not shown, 
for it is much smaller than the estimated contact term 
contribution. In addition, we argued in Sec. 13.11 that we 
do not expect large corrections to the ir-^He scattering 
length due to the net effect of the dispersive and the A 
contributions. Thus, our prediction for the 7r-^He scatter- 
ing length is 

^(m+2iv+3Ar) ^ ± 4 ± 7) X m-'m-' (39) 

where the first uncertainty is due to ambiguities in tt- 
N scattering lengths whereas the second one represents 
the uncertainty in the few-nucleon effect. This result does 
not include isospin violating few-nucleon effects. Those are 
found to be sizable for tt-^H scattering [77] although signif- 
icantly smaller than the theoretical uncertainty of 7r-^He 
calculation. The result of Ea. ((M)) is to be compared to the 

Table 8. Experimental results for the vr-'^He scattering length. 
The entries contain measured energy level shifts of vr-'^He 
atomic bound states together with corresponding scattering 
lengths. 

eis [eV] a^_3He [lO~^ m"'] 

R. Abela et al. ^ 44 ± 5 56 ± 6 
G. R. Mason et aZ. [79] 34 ± 4 43 ± 5 
I. Schwanner et al. [gO] 32 ± 3 41 ± 4 



experimental results for the tt-'^Hc scattering lengths given 
in Table [5] Those are extracted from the measurements of 
the Is level shifts in 7r~-'^He atom due to the strong inter- 
actions [7511751 150] by using DGBT-type formulae [5T] and 
including logarithmic corrections of Ref. [55]. The table 



Table 9. Numerical results for the relative scalings of the few- 
nucleon contributions. 





A/ A 


",r-3Ho/"7r-3Hc 


(3N) , (2N) 


CD-Bonn 
AV18 




0.220 
0.222 


0.221 
0.194 


NLO 
NLO 
NLO 
NLO 
NLO 


400/500 
550/500 
550/600 
400/700 
550/700 


0.206 
0.212 
0.210 
0.205 
0.210 


0.227 
0.181 
0.175 
0.228 
0.180 


N^LO 
N^LO 
N^LO 
N^LO 
N^LO 


450/500 
600/500 
550/600 
450/700 
600/700 


0.212 
0.220 
0.215 
0.207 
0.216 


0.222 
0.194 
0.224 
0.236 
0.224 



demonstrates that the results of the first measurement are 
in contradiction with the others even within the large ex- 
perimental uncertainties. It is getting even more intriguing 
because it is only the first measurement that agrees with 
our theoretical prediction ([M]) based on ChPT. Clearly, 
a new measurement of this quantity recently performed 
at PSI [5^ is of high importance to resolve the existing 
discrepancies. 

In Table [HI we have compiled the relative contribu- 
tions of one-nucleon, two-nucleon, and three-nucleon dia- 
grams (a^^i^ijj^, a|f^4^, and a'^^^^^J. We omitted the re- 
sults of the LO wave functions here, since their description 
of the binding energies is generally poor. It sticks out that, 
numerically, the suppression of few-nucleon corrections is 
less than expected by Weinberg's power counting. Based 
on these findings, we are led to the conclusion that the 
power counting gives reasonable guidance in identifying 
the most important contributions, however, for a quan- 
titative understanding, explicit calculations for the lead- 
ing few-nucleon contributions are necessary to estimate 
the contribution of the class of iV-nucleon contributions. 
Qualitatively, more-nucleon diagrams are still sufficiently 
suppressed so that the series of one- , two-, ... nucleon con- 
tributions can be truncated at sufficiently low complexity 
of the problem — we find a factor of 5 suppression when 
going from an iV-nucleon operator to an {N + l)-nucleon 
operator. It is important to note in this context that we 
find the same suppression factor for = 1, = 2, and, as 
will be shown in the next section, A^ = 3. Due to this, the 
four- nucleon diagrams turn out to be already insignificant. 

In summary, we have calculated Tr-'^He scattering length 
including leading three-nucleon terms and two-nucleon op- 
erators. Due to the presence of the large contact term 
in the isovector channel the present calculation basically 
reaches the edge of the theoretical accuracy. We also find 
that the results are in qualitative agreement with Wein- 
berg's counting. 
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Table 10. Summary of the shifts of the 7r-^He scattering length due to the few-nucleon corrections '^^_4fjp) ^tt— ^He' ^tt— ^He' 
and a^4jj^. For NLO and N^LO, both cutoffs {A/ A) are given in MeV. Central values and standard deviation for the scattering 
length results are given in lQ~^m^^ . 





A/ A 




„|I^^) 






AV18 




-49.5(7) 


-1.29(2) 


5.00(5) 


2.73(84) 


NLO 


400/500 


-56.1(15) 


3.02(1) 


7.18(1) 


2.79(11) 


NLO 


550/500 


-51.0(8) 


0.41(1) 


6.04(2) 


2.18(22) 


NLO 


550/600 


-51.4(5) 


0.53(1) 


5.95(1) 


2.02(22) 


NLO 


400/700 


-54.5(3) 


3.10(1) 


7.13(2) 


3.92(42) 


NLO 


550/700 


-51.6(12) 


0.58(1) 


5.72(1) 


2.54(59) 


N^LO 


450/500 


-54.4(4) 


1.92(1) 


6.98(2) 


3.00(20) 


N^LO 


600/500 


-52.0(8) 


-0.09(2) 


6.16(3) 


2.13(11) 


N^LO 


550/600 


-52.7(6) 


0.50(1) 


6.42(3) 


2.31(31) 


N^LO 


450/700 


-52.7(7) 


1.81(1) 


6.68(3) 


2.56(10) 


N^LO 


600/700 


-53.9(8) 


0.36(1) 


6.34(2) 


2.81(17) 



8 Two- and four-nucleon contributions to 

Finally, we want to discuss the few-nucleon contributions 
to 7r-'*He scattering. Because of the isovector character of 
the leading three-nucleon contributions, we here only need 
to study two- and four-nucleon ones. Our results are sum- 
marized in Table [TOl Qualitatively, the results for the two- 
nucleon operators are similar to the ones for tt-^H and for 

7r-"^He scattering. The leading two-nucleon term a^ljjp is 
the by far most important contribution depending on a*-"-* . 
Just based on the symmetry arguments al^ljj^, should be 
approximately twice as large as the one in tt-'^Hc scatter- 
ing. As one can see from Tables [7] and [TUl this is indeed 
supported by our results. In addition, the contribution 
a^ljjg is approximately twice as large as for the deuteron 
although the number of NN pairs is six. This can be ex- 
pected from the isospin structure of this amplitude leading 
to opposite signs for the contributions of neutron-proton 
pairs and proton-proton (or neutron- neutron) pairs [84,. 
Trivially, the one-nucleon terms in tt-^Hc and in tt-^H scat- 
tering also scale with a factor 2, as seen from Eq. (P). 
Unfortunately, this implies that the correlation of a'^^^ 
and a'^^^ due to experimental results for 7r-^He scatter- 
ing is very similar to the one based on a tt-^H analysis. 
The most important correction to the leading two-nucleon 
term originates from the triple scattering diagram. Again, 
it is enhanced compared to the estimate of naive dimen- 
sional analysis. In full analogy with 7r-^He and tt-'^Hc scat- 
tering the ratio a^4fjc/'^^4fjc about 0.1-0.13 and is 
only smoothly dependent on the cutoff. 

From the cutoff variation of the N^LO results, we es- 
timate that missing counter terms should contribute of 
the order of 2 x lQ~^m~^. This is approximately 4 % 
and in line with our expectations from Weinberg count- 
ing. a^4He again suppressed and negligible compared 
to the cutoff dependence of a^ljj^. Similar to the results 
for TT-^He, even the sign of this contribution is not fixed 



clearly indicating the sensitivity to the short distance part 
of the wave functions. 



Let us now turn to the contribution of the four-nucleon 
operator aj^4jjg. Its contribution is opposite in sign to the 
leading two-nucleon contribution, but similar in size as the 

cutoff variation of a^lfj^, and, therefore, comparable to 
the contribution of the first isoscalar two-nucleon counter 
terms. The relative suppression of the a^4jjg compared 

to a^ljjg is approximately 0.06. From Weinberg counting, 
however, we naively expect a much larger suppression of 
(m.n-/mAr)^ ~ 4 X lO"**. On the other hand, this deviation 
was expected and in line with the results of our calcula- 
tions for TT-^He where we also observed a significant overes- 

timation of the calculated ratio o^sho/'^^shc compared 
to the dimensional analysis. In fact we find numerically 
a suppression of about 1/5^ for the four-nucleon operator 
compared to the two-nucleon operator, completely in line 
with the pattern of successive suppression found in case 
of TT-^He scattering. As a result of this, the contribution 
of the four-nucleon operators can be, in principle, omit- 
ted, for it coincides in the magnitude with the estimated 
counter term contribution. This is the reason why we do 
not take into account the other four-nucleon topologies 
analogous to those given in Fig. [T3] for tt-^Hc scattering: 
Based on the study of the tt-^Hc process, we expect their 
contribution to be suppressed as compared to the leading 
four-nucleon operator. 



Thus, we conclude that pion scattering on "^He can 
not provide any additional constraints on a^+^ and a^"-* 
although this process will be helpful to study the system- 
atic uncertainties of experiment and theoretical calcula- 
tions. Again we find that Weinberg counting is not able 
to capture quantitatively the relative suppression of oper- 
ators with a different number of active nucleons. 
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9 Conclusions and Outlook 

In this paper, we have systematically studied the few- 
nucleon contributions to pion-nucleus scattering. We com- 
pared our numerical findings to power counting estimates 
of two power-counting schemes, Weinberg counting and 
Q-counting. Our numerical results do qualitatively sup- 
port Weinberg's power counting but not Q-counting. An 
evidence for this was provided by studying the dependence 
of different contributions on the binding energy of the 
deuteron. This dependence is expected to be significantly 
different in these two counting schemes varying from a 
constant (up to logarithmic corrections) for the ratio of 
relevant two-nucleon contributions in Weinberg counting 
to a power law behavior in Q-counting. It is shown both 
numerically and analytically that the ratio is indeed very 
weakly depending on the binding energy in clear favour 
of Weinberg counting. In the course of the analytic analy- 
sis, we found that wave functions based on a perturbative 
treatment of pious should not be used for the calculation 
of pion-nucleus scattering. Especially, a treatment of the 
pion as a heavy field would be justified only for unphysi- 
cally small deuteron binding energies below 0.1 MeV. 

Numerical results for some particular diagrams are not 
in line with the estimates based on Weinberg counting. 
This was resolved when we could trace back enhancement 
or suppression of the pertinent diagrams either to acciden- 
tal spin-isospin factors or to specific properties of the loop 
function of the triple scattering diagram, which can not 
be, of course, taken into account in the counting scheme. 
Especially, the last insight might be of relevance for the 
treatment of pion production and also NN forces. 

Generally, our results point to a smaller suppression 
of more-nucleon terms than expected by naive dimen- 
sional analysis. To be specific, our results demonstrate 
that the ratio of one-nucleon to two-nucleon to three- 
nucleon to four-nucleon contributions scales roughly as 
100:25:5:1 as compared to the pattern 50^:50^:50:1 which 
can be expected from Weinberg counting. Apparently, the 
additional integral measures that enter when nucleons are 
added and that provide most of the suppression in the di- 
mensional analysis are partially canceled by other mecha- 
nisms. It should be stressed that even the relatively mild 
suppression found is still quantitatively sufficient to allow 
for controlled calculations for pion-nucleus scattering. Es- 
pecially, for the TT-^He scattering length the four-nucleon 
operator is already numerically insignificant. 

We performed explicit calculations up to four-nucleon 
contributions to pion-nucleus scattering. Based on these 
calculations, we conclude that three-nucleon contributions 
might still be relevant. On the other hand, it turned out 
that the most important three-nucleon diagram is isovec- 
tor. Therefore, for the important case of ^He, only one- 
nucleon and two-nucleon diagrams will contribute signifi- 
cantly. It turned out that the relative contribution of one- 
and two-nucleon diagrams is very similar for tt-^H and tt- 
^He. Therefore we expect that the corresponding bands 
a'-"'"'' vs a*^") for tt-^H and 7r-^He will be almost on top 
of each other. Thus, we tend to conclude that pion scat- 
tering on '*IIe can not provide any additional constraints 



on a^"*"^ and a''~\ however, an improved measurement of 
the 7r~-^He atom will provide an additional cross check 
of the systematics of the analysis. On the other hand, tt- 
■^He does contain a new non-trivial dependence on a'"-'. 
That is why we have studied this process in more de- 
tail. For the first time three-nucleon contributions were 
included. Based on the latest numbers for a'^-* — which 
denotes the isoscalar scattering length including the lead- 
ing, isoscalar isospin violating corrections — and a^^^ we 
obtained a^_3Ho = (62±4±7) x 10~'^m~^ where the uncer- 
tainties are due to ambiguities in tt-N scattering lengths 
and in the unknown isovector two-nucleon contact opera- 
tor, respectively. Unfortunately, the large theoretical un- 
certainty related to the contact term in the isovector chan- 
nel precludes one from considering Tr-'^He scattering on the 
same footing with tt-H and tt-^H processes to extract a^~^^ 
and a*-"-* . On the other hand, the results of a new measure- 
ment of TT-'^He at PSI ^5] which are currently analyzed, 
are of high interest to check the predictions of ChPT, espe- 
cially since the theoretical value lies above the most recent 
experimental value (we find a 2a discrepancy), while it is 
consistent with an earlier measurement. 

Our results also showed that we are able to predict 
few-nucleon contributions for pion scattering off isoscalar 
(isovector) nuclei up to a level of 5 % (30 /^El) of the leading 
two-nucleon contributions before additional counter terms 
enter. Since isospin violation is predicted to contribute at 
the few-percent level as well, it will be required to further 
work on IV contributions, before we finally extract a^"*"^ 
and a^~^ from the available data on pionic atoms. Work 
in this direction is in progress [77] . 

We thank E. Epelbaum, D. Gotta, M. Hoferichter, B. Kubis, A. 
Kudryavtsev, U.-G. Meifiner, and D. R. Phillips for useful dis- 
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HadronPhysics2 project. The work of V.B. was supported by 
the State Corporation of the Russian Federation "Rosatom". 
Part of the numerical calculations have been performed on the 
super computer cluster of the NIC, Jiilich, Germany. 

A Independence of three-nucleon diagrams 
on the parameterization of the pion field 

In this appendix, we explicitly show the independence of 
the sum of Feynman diagrams depicted in Fig. [T3]on the 
parametrization of the pion field. For this, we closely fol- 
low the recipe outlined in [S^- We start with the most 
general expression for the chiral matrix U (tt) 

C/(7r)=exp(^^(r.7r),g(7rV./^)). (40) 

where the function g is an arbitrary regular function with 
the property g (0) = 1. After expansion up to second order 

This is 10% of the full scattering length, since for isovector 
targets, contrary to their isoscalar counter terms, the scattering 
length is dominated by the one-nucleon term. 
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in the pion fields it can be written as 



1 

6 



f2 



(41) 



For a = — i this equals the expression for the chiral matrix 
U in the so-called exponential gauge, U — exp (^j-iT- tt)^ 

-T ■ TT, which is 



In the cr-gauge we have U = — j2 i j 

reproduced up to terms of fourth order in the pion fields 
by using a = (for details see [5^H5] ). 

Based on this ansatz for the chiral matrix U, it is easy 
to derive the corresponding Feynman rules for the NNtt, 
NN27r, NNStt, NN47r, and An vertices. The first two of the 
vertices, turn out to be independent of a: 



VnNti 
VnN2tt 



9A 

1 

4/2 



(T-9)r'^ 



2«- (92 



qi)e'''"=T'' 



(42) 
(43) 



The third, fourth, and fifth one depend on a. They are 
given by 



9A 



{^'^V'^ [qi + 92 + 4a (qi + 92 + Qs)] 
[qi + qs + Aa {qi + 92 + Qs)] 
hS'^'r" [q2 + 93 + 4a {qi + 92 + 93)] } (44) 



V, 



1 

«74 



Ba) V 



•{^'^V^'^ (94-93) + V'^'' (94 
+^^V'^''(g2-'Zi) + '5'^V"^ (94 



91) 
92) 



_gad^fbc ^ 




For these expressions, we assume that all momenta are 
going out of the vertices. 

Here, we are not interested in the expressions for the 
amplitudes, but only in the a dependence. Therefore, for 
simplicity, we only keep the a-dependent terms and show 
that these vanish for the sum of the diagrams in Fig. [141 

We identified seven a-dependent diagrams in leading 
order, which we summarize in Fig. 1141 Choosing v — 



(1,0,0,0), it is easy to derive for the a-dependent part 
of the threshold amplitude shown in Diagram (a) 



iM 



Ell) 



-a 



4/6 {qf - ml) {ql - m^) 
(0-1 • qi) (0-3 • qr3)r2^ 



I gad/ (^^b^d ^ ^d^b^^ ^ ^dbf (^^a^d 



+26'^''^ (n • Ts)} 

Similarly, for Diagram (b)-(g) one finds 



iM^^^a^-^{rT,-q,){a,-q,)ri 

1 1 



(47) 



+ 



(9? - ml) {ql - ml) {qf - ml) {ql - ml) 
^ ^ (48) 



m - ml) {qi - ml) 



2 

a (o-i • qi) (0-3 • 93) T, 



4/.6 

{^adf ^^b^d^^d^b^^^abf (^r,-Ts)} 



1 



1 



iM 



Gil) 



(9? - ml) [ql - ml) (qf - ml) (gf - ml) 
{qi-ml){ql-ml)\ ^"^^^ 



dAmn 



4/6 {qf ~ ml) (qj ~ ml) 
(cTi • qi) (0-3 • qs) tI 

X [e'^'^f {t^tI + rtr^) + e^'f {t, ■ t,) }(50) 



iM 



9\m^ 



4/6 {ql-ml){ql-m:i) 
{a I ■ qi) (0-3 • 93) T2^ 

,ad/ (hd , _rf_b^ I ,ab/ 



{e'^^f {rlr^ + r^rl) + e^'^ [r, ■ r,) }(51) 



g\m^ 



4/6 (ql _ ^2 ) (^2 _ ^2 ) 



(cTi • qi) (0-3 • 93) r2^ 

X [e'^^'f (tM + T^T^) + e'^^'f (n • T3) }(52) 



g\m^ 



4/6 (92 _ „4) (g2 _ ^2 ) 
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(a) 



— • — 

c T qi 



/ 

'2 

d i qs 



(b) 



b ,-'^Ti 



Q2 



d i 



Q3 



<72 



T3 



T3 



(d) 



Tl 



(f) 



Tl 



\ 92 



\ / q2 



92 \ 



93 



T3 



T3 







• 




a 


^^^W i 






^2 ; d 






• 



T3 
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Fig. 14. Individually a-dependent leading three-nucleon contributions to pion-nucleus scattering. 



{(Tl ■ qi) ((T3 • qj,)Tl 



Tl Tq 



i'3 T- '1 '3; -i-e"-^ i-Ti •T3)}(53) 



Based on these expressions, it is easy to convince oneself 
that the sum of the expressions Eqs. (|T7)) to ([55)1 cancels. 
This implies that the sum of the diagrams is independent 
of a. 



B Explicit expressions for the 
half-Coulombian three-nucleon diagrams 

In this appendix, we summarize the explicit expressions 
for the amplitudes of Fig. [141 The sum of these diagrams 
is independent of the parametrization of the pion field as 
discussed in Appendix [X] We refer to these amplitudes as 
half-Coulombian because for most of the diagrams of this 
group one of the pion propagators is . They read 



32/6 (q2+^2)(q| 



X {(Tl ■ Qi) (0-3 • 93) 



I T, 



+2e'''f (ti-T3)}t; 



/ 
2 



iM 



Git) 



Qi iQ?+ml) {qi 



X (cTi • 9i) (era ■ 93) 



x{e""-' (qf TfT- 



2 ^ad 
1 T^X^-i 



iM 



llMl ') 



(92 



8/6 q|(q2 + ^2)(q2+^2) 

X (cri ■ 91 ) (cr3 ■ 93) 

x{e'^'f{q^ir'iri + qirtrl) 
+ {qi + ml)e'^'f{Ti-T,))ri 
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iM 



1 



(cTi • gi) 



16/6 qi{qf+ml) 



3 J 



iM 



92 e"^"^ (ri • T3)}t, 



f 



16/^ qiiqi + ml) 



iM 



Git) 



+92 e""^ (ri • r3)}r, 



f 



+92 e^'f (xi • T3)}t/ 



(57) 



16/6 qi{qf + ml) 
+92 e'^'f (n • X3)}r/ (58) 



(59) 



(60) 

The definition of the momenta can again be read off from 
the figures. For these expressions, we have assumed a = 0, 
which corresponds to the a-gauge. 



where q = \/ p"^ + p'"^ — 2pp'x is an internal momentum 
transfer and Pk{x) is the degree k Legendre polynomial. 

With these definitions the PW decomposition of the 
amplitudes Eqs. dH), Q and dH]) read 



{a'\M 



i-y 



(64) 



^2 ^2 



1 1 



a'|A^«^=^|a) = |5|E E (2/ + DM no 



//o 



(11/, 000), 



p'^\~pr -giip'p) 



(2Aii)!(2^2 



k k 0\ 

{2k + 1) <( Ml M2 / > {k^ill', 000) {k^i2l 000) (-)' 

^' ^ / J 

s' s / ^ ( V I f 



18v/(2s + l)(2s' + l)(2j + l) Is' sf 



i i 1 i' j 



(65) 



C PW decomposition of the two-nucleon 
operators 

In this appendix, we briefly summarize the PW decom- 
posed expressions for the pertinent two-nucleon operators. 
The PW states for the NN system read 

I") = \p{ls)jm) (61) 

where p is the magnitude of the NN relative momentum, 
I and s are the corresponding orbital angular momentum 
and NN spin and j,m are the total angular momentum and 
its projection. We abbreviate this set of quantum numbers 
by a. Here we only consider tt^-^H scattering. Therefore, 
the two nucleons are in an |t = Omt = 0) isospin state 
and the cartesian components a,b = 1...3 of the pion are 
(l/\/2,-i/\/2,0). The pertinent isospin matrix elements 
are 

{t ^ Omt ^ 0|(5°'' (ri • T2) \t = Omt = 0) = -3 
{t = Omt = 0\t'It^ + T'^T^lt = Omt = 0) = -2 

(62) 

The spin-orbital part can be expressed in terms of in- 
tegrals 

gUpp)^M-)''V2k + T I'^dxPkix)^ 

~giip'p)^2ni-rV2kTT j\x Pk{x)^j-^^-^ 

(63) 



{a'\M^\a} 

(66) 

for vr^-^H scattering in terms of Qj-coefficients and Clebsch- 
Gordan coefficients (J1J2J3, 771177127713). Here, we normalize 
the momentum eigenstates such that the expectation val- 
ues read 

{M) 

poo 

= EE / dpp'^ dp' p'^ 'ipi,^,y{p'){a'\M\a)ipisj{p) 

(67) 

for wave functions normalized to 
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